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o ; Abstract 

The activated rate process for non-equilibrium open systems is studied taking into account both 
internal and external noise fluctuations in a unified way. The probability of a particle diffusing 
,-SZh " passing over the saddle point and the rate constant together with the effective transmission coeffi- 

cient are calculated via the method of reactive flux. We find that the complexity of internal noise 
is always harmful to the diffusion of particles. However the external modulation may be beneficial 
to the rate process. 
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I. INTRODUCTION 



More than seventy years ago H. A. Kramers published a seminal work on the diffusion 
model of chemical reactions Ever since then, the theory of activated processes has become 
a central issue in many fields of study 

aa 

, notably in chemical physics, nonlinear optics 
and condensed matter physics. In the model, the particle was supposed to be immersed in 
a huge equilibrium medium so that it can gain enough energy to cross the barrier from its 
thermally activation. The common feature of an overwhelming majority of such treatments 
is that the system is thermodynamically closed. That is to say, the noise of the medium 
is of internal origin so that the dissipation-fluctuation theorem 0, 5} is satisfied and a zero 
current steady state situation is characterized by an equilibrium Boltzmann distribution. 

However, when the system is thermodynamically open, for example, driven by an external 
noise which is independent of the medium [6|, no relation between the dissipation and 
fluctuations can be dependent on. The corresponding situation aforesaid, if attainable, can 
then only be defined by a steady state condition 3, S] which may depend not only on the 
strength and correlation of the external noise but also on the dissipation of the system. The 
external noise will modify the dynamics of activation in the region around the barrier top 
so that an unusual effective stationary flux across it gets resulted. This would in no doubt 
induce an unfamiliar activated barrier escaping process that is worth pursuing for the rate 
theory. 

Therefore we present in this paper a recent study of us on the activated rate process 
where the diffusing particle is under the joint influence of the internal noise combining an 
external one. The paper is organized as follows: In Sec. II, reactive dynamics at the barrier 
top is investigated by analytically solving the generalized Langevin equation. In Sec. Ill, we 
give a detailed discussion about the combined effect of internal and external noises on the 
rate process by asymptotically calculating the rate function and its transmission coefficient. 
Sec. IV serves as a summary of our conclusion where some implicate applications of this 
study are also discussed. 
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II. REACTIVE DYNAMICS AT THE BARRIER TOP 



We consider the motion of a particle of unit mass moving in a Kramers type potential 
U(x) such that it is acted upon by random forces ((t) and e(t) of both internal and external 
origin, respectively, in terms of the following generalized Langevin equation (GLE): 



x + / dt , 7(t-t')x(t') + d x U(x) = ((t) + e(t), 



(1) 



where U(x) is the potential, the friction kernel l{t) is connected to internal noise by the 



well-known fluctuation-dissipation theorem (FDT) [4], |5[ (((t)((t')} = kBT^(t — £'). Both 
the noises ((t) and e(t) are assumed stationary and Gaussian with arbitrary decaying type 
of correlation. We further assume, without any loss of generality, that ((t) is independent 
of e(t) so that we have (C(t)e(t')) = and (e(t)e(t')) e = 2Dip(t - 1'). Here (• • • ) e implies the 
averaging over all the realizations of e(t) with D the intensity constant and ip(t) a relevant 
memory function. In other words, the external noise is independent of the friction kernel 
7(t) and so there is no corresponding fluctuation-dissipation relation. However, correlation 
{e(t)e(t')) e = 2Dip(t — t') is reminiscent of the familiar FDT formula due to the appearance 
of the external noise intensity, it serves rather as a thermodynamic consistency condition 
instead. 

Due to the Gaussian property of the noises ((t) and e(t) and the linearity of the GLE, the 
joint probability density function of the system oscillator must still be written in a Gaussian 
form [9] The reduced distribution function can then be yielded by integrating out all the 
variants except x as 



W(x,t; x ,v ) 



1 



-exp 



(t) 



(2) 



V27T a x (t) 

in which the average position (x(t)) and variance cx^(t) can be obtained by Laplace solving 
the GLE. In the case of an inverse harmonic potential U(x) = —^mulx 2 , it reads 



(*(*)) 



1 + ^ 2 



H(t')dt' 



dt x H{t - ti 



x + H(t)v 
dh{i{t x )t{t 2 ))H{t-t 2 ) 



(3a) 
(3b) 



where *(*) namely the response function can be vielded from inverse Laplace transforming 



His) 



_1 with residue theorem [lCH. Illfl. t;(t) (= ((t) + e(t)) is an effective 
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noise of zero mean (£(£)) = whose correlation is given by 

mat')) = mat')) + <e(f) C (o>., (4) 

where the two averages in the right hand side are taken independently. 

The probability of passing over the saddle point, namely also the characteristic function 
which is crucial for the activated barrier crossing process, can then be determined mathe- 
matically by integrating Eq.([2J) over x from zero to infinity as 



P(x 0l v ;t) = I W(x,t;x 0} v )dx, 

(x(t)) 







-erfc 
2 



V2a x (t) ' (5) 
The escape rate of a particle, defined in the spirit of reactive flux method by assuming the 
initial conditions to be at the top of the barrier, can then be yielded from 
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— OO 



k(t) = t I dx o I v W st (x , v Q )P{x , v ; t)5(x - x b )dv (6) 



in the phase space. This in proceeding results in a generalized transition state (TST) rate 
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14|: k TST = u b/( D b+^(°°)) an d an effective transmission factor 



2(.i\ \ - 1 / 2 

«(*)=(! + > (7) 



D b H 2 (t) 

where W s %{xq,v ) = ^exp[— {^- + d~^t^)}} is a Boltzmann form stationary probability 
distribution which can be obtained from the steady state Fokker-Planck equation jg]. This 
stationary distribution for the non-equilibrium open system is not an equilibrium distribution 
but it plays the role of an equilibrium distribution of the closed system, which may, however, 
be recovered in the absence of the external noise. (D& + \&(oo))/ kg in the exponential factor 
of k TST defines a new effective temperature characteristic of the steady state of the non- 
equilibrium open system. 

In the particular case we have considered, parameters used heretofore are defined as [?], [s| : 
U(x) = Ui, — \Q 2 {x — Xf,) 2 the renormalized linear potential near the barrier top with Q b 
an effective frequency and £/& the barrier height, ^(oo) and D b = <3>(oo)/T(oo) are to be 
calculated from 

$(*) = nl(t)a 2 xv (t) + T(t)a 2 (t) + (8a) 
9{t) = j t a 2 xv {t) + nt)°l(t) + n 2 (t)a 2 x (t) - a 2 v (t), (8b) 



for the steady state in which T(t) = -|lnA(t), tlftt) = H(t)(H(t) - H(t))/A(t) and A(t) = 
^-[l-^£drH(r)] + H 2 (t). Other variances besides cr^(t) are also to be got from Laplace 
solving the GLE. These variables will play a decisive role in the calculation of barrier escaping 
rate. Therefore, in general, one has to work out these quantities first for analytically tractable 
models js]. 

As is expected, all the parameters besides the rates n(t) and k(t) are closely related to the 
internal and external noises. Therefore a combining control of internal and external noise 
on the activated barrier escaping process is prospected. This is what will be involved in the 
following sections. 

III. INTERNAL VS EXTERNAL NOISE 

Before accomplishing the following calculations, let us firstly digress a little bit about 
P(xq, vq; t) and n(t) which are the central results of this study. As has been shown in Eqs. (J5j) 
and (171), both the expressions of P(xq, v ; t) and n(t) are reminiscent of the familiar previous 
results [l^Q]- Although variance cr^(t) has been changed intrinsically by the external noise, 
difference lives only superficially in the emergence of which is an asymptotical constant in 
the long time limit. Due to the independence of internal and external noises, other variances 
such as (x(t)) and H(t) depend only on the internal noise. Therefore, from the viewpoint 
of diffusing passing over the saddle point, the mean position of the Gaussian packet relies 
simply on the internal noise while the width of it depends not only on the internal noise 
but also on the external one. It is the combining effect of internal and external noise that 
determines the final diffusing process. In what follows, it shall be concerned with several 
limiting situations to illustrate the general result systematically for both thermal and non- 
thermal activated processes. 

A. Internal white noise 

Firstly we consider the simplest case of a 5 correlated internal thermal noise combining 
with no external ones. To this end, it is to set 

e(t) = and (C(t)C(O) = k B T^8(t - t'). (9) 
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By combining with the abbreviations in Eqs. (l3~]) and ©, all the quantities we need for the 
activated process can be obtained easily. After some algebra it follows that 

*(oo) = and D b = k B T. (10) 

The relations obtained heretofore reduce to the general form 

^-(^■AtY*. (id 



k B TH\t), 

This is a trivial result for the one- dimensional time-dependent barrier passage [16|. It 
generally describes the possibility of a particle already escaped from the metastable well to 
recross the barrier. In the case of no external noise modulation, it keeps its usual form just 
as it should be. 

B. Internal color noise 

Next we discuss the case of purely internal color noise. For example, we set the internal 
noise to be Ornstein-Uhlenbeck (OU) type Q. That is 

e(t) = and (C(t)C(O) = A^T-e^*' 1 ^, (12) 

where G denotes the strength while r c refers to the correlation time of the noise. It should 
be noted that for r c — > the internal noise shown above becomes also 5 correlated. After 
some algebra we can obtain from Eqs.Q and © again that 

tf(oo) = k B T ( % - l) while D b = k B T. (13) 



This will result in a different form of TST rate k TST but has no influence on the form of 
transmission coefficient n(t). However, since in most cases = uj%, the value of \l/(oo) 
is actually close to 0. The values of variance cr^(t) and H(t) are also changed intrinsi- 
cally. Therefore, although the form of n(t) has not been changed, the rate process has been 
modified implicitly due to the alteration of internal noise. 

C. Internal and external white noise 

In further, let us turn to investigate the more complicated combining case where both 
the internal and external noise to be 5 correlated, i.e. 

(e(t)e(t')) e = 2V5(t - t') and (C(*)C(*')> = k B Tj5(t - t'). (14) 



in which T> is the strength of the external white noise. This may be the simplest combining 
case of internal and external noise. Derive again from the abbreviations heretofore we obtain 

V 

tf(oo)=0 and D b = k B T + — . (15) 
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Noticing that comparing with the aforesaid case in Sec JIII Al a new effective temperature due 
to external noise is defined here by 2?/ (7/15) [19]. But in the limit of T> — Y 0, the previous 
one- dimensional form of n(t) for pure internal white noise can still be recovered by Eq.(J7j). 
However, the rate process has also been changed intrinsically due to the effect of external 
noise. 

D. Internal color and external white noise 

Finally, we consider a particular case where the external noise is S correlated while the 
internal is an OU process, i.e., 

(e(t)e(t')) e = 2V5(t - f) and (C(t)C(O) = A^T-e^*'^, (16) 

both symmetric with respect to the time argument and assumed to be uncorrelated with 
each other. By virtue of similar derivations as heretofore, we find it is difficult to get an 
explicitly simple expression of ^(00) while Df, = k B T is recovered again. This is nontrivial 
in the presence of external noise because it seems as if there is not any effect on the rate 
process that comes from the external noise. But actually all the effects have been contained 
in the calculations of cr^(t)) and H(t) so as in that of n(t). 

In order to give an explicit revelation of the combining effect of the noises on the rate 
process, we plot in FigJTJthe instantaneous values of P(x ,v ;t) and n(t) at different types 
of combining cases aforesaid. From which we can see that the asymptotic stationary value of 
P(xq, Vo] t) and K,(t) (defined as P st and K st respectively) in the internal color case is smaller 
than that of internal white case no matter the system is modulated by an external white noise 
or not. On the contrary, both P st and n st are larger than those of pure internal case supposing 
an external noise is set on with modulation. Thus we can infer from considering the intrinsic 
meaning of P(xq, Vq] t) and f«(i) that the complexity of internal noise (or dissipation) is always 
harmful to the diffusion of particles. However the external modulation may be beneficial 
to the rate process. This is a non-trivial result of great meaning to many different kinds of 
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FIG. 1: Instantaneous values of x( x 0, v 0'it) an d K>(t) for various types of internal and external 
noise combinations. Dimensionless parameters used are 7 = 4.0, uib = 1.0, G = 3.0 and r = 0.3. 
Particles are assumed to start from position xq = —1.0 with initial velocity v$ = 6.0. Lines signed 
as A, B, C, D represent the different cases discussed in the corresponding subsections respectively. 



realistic situations in forming a non-equilibrium (or non-thermal) system-reservoir coupling 
environment 
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IV. SUMMARY AND DISCUSSION 



In summary, we have studied in this paper the activated rate process for non-equilibrium 
open systems taking into account both internal and external noise fluctuations in a unified 
way. We calculated the probability of a particle diffusing passing over the saddle point 
and the rate constant together with the effective transmission coefficient via the method of 
reactive flux. The combining control of internal and external noises on the activated barrier 
escaping process is investigated. We find that the complexity of internal noise is always 
harmful to the diffusion of particles. However the external modulation may be beneficial to 
the rate process. 

We believe that these considerations are likely to be important in other related issues 
in non-equilibrium open systems and may serve as a basis for studying processes occurring 



within irreversibly driven environments 22|, |23( and for thermal ratchet problems 24j . The 



externally generated non-equilibrium fluctuations can bias the Brownian motion of a particle 
in an anisotropic medium and may also be used for designing molecular motors and pumps. 
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